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Strong rate of convergence for the Euler-Maruyama 
approximation of SDEs with Holder continuous drift 

coefficient 

Olivier Menoukeu Pamen* and Dai Taguchd 


Abstract 

In this paper, we consider a numerical approximation of the stochastic differential equation 
(SDE) 



where the drift coefficient b : [0, T] x R d —> R d is Holder continuous in both time and space 
variables and the noise L = ( L t )o<t<r is a d-dimensional Levy process. We provide the rate 
of convergence for the Euler-Maruyama approximation when L is a Wiener process or a trun¬ 
cated symmetric a-stable process with a G (1,2). Our technique is based on the regularity of 
the solution to the associated Kolmogorov equation. 
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1 Introduction 


Let X = (X t )o <t<T be the unique strong solution to the following d-dimensional SDE 


Xt — xo + f b(s , V s )ds + Lt, Xq £ t G [0,T], (1) 

Jo 

where L = (L t )o<t<T is a d-dimensional Levy process on a complete probability space (fi, T, P). 
The drift coefficient b : [0,T] x —> W 1 is assumed to be 77 -Holder continuous in time with 

77 G [1/2, 1] and bounded and /3-Holder continuous in space, i.e., 


IhllcVfOTl) := SU P \ b (t,x)\ 

b te[0,T],a:eR d 


I b(t,x) - b(t,y) I 

te[o S r] P ^y \x~y\ p 


< OO. 
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Here /3 satisfies some conditions (see Theorem 12.111 and 12.131) . 

Consider the Euler-Maruyama approximation of SDE JT]) given by 

X t (n) = Xo + J o b (th(s), ^j t) ) ds + L t , t G [0, T], (2) 

where r/ n (s) = kT/n =: t^ if s € [kT/n, (k + l)T/n). It is well-known (see for example fill ) 
that if L is a Wiener process and the coefficient b is Lipschitz continuous in space and 1/2-Holder 
continuous in time then the Euler-Maruyama scheme has strong rate of convergence 1/2, i.e. for 
any p > 0, there exists C p > 0 such that 


sup 

x t - xj n) 

P~ 

.0 <t<T 


_ 


< 


C p 
n p / 2 


In general, convergence of numerical methods are studied under the above assumption that 
the coefficient of the SDE is globally Lipschitz continuous. However, this assumption is not always 
satisfied for SDEs used in practice (for example mathematical finance, optimal control problem 
and filtering), this makes the global Lipschitz based idea not immediately applicable. This lead to 
the study of Euler-Maruyama approximation for SDEs with irregular coefficients. 

When L is a Wiener process and & is a continuous function satisfying a linear growth condition, 
Kaneko and Nakao [ 8 ] proved that the unique strong solution of the SDE (JT]) if it exists can be 
constructed as the limit of the Euler-Maruyama approximation. However, they do not study the 
rate of convergence for the Euler-Maruyama scheme. Hashimoto [7] considers Euler-Maruyama 
approximation for solutions of one dimensional SDE driven by a symmetric a-stable process. Under 
the Komatsu condition (a Holder type condition), the author shows strong convergence of the 
Euler-Maruyama approximation. The proof of the result is based on an approximation argument 
introduce in [T2| and which is very close to the Yamada-Watanabe approximation technique. 

In the recent years, there has been a lot of studies on the rate of convergence for the Euler- 
Maruyama approximation with non-Lipschitz coefficients. In the diffusion case, assuming that 
the drift coefficient is the sum of a Lipschitz continuous function and monotone decreasing Holder 
continuous function, [5] gives the order of the strong rate of convergence for one-dimensional SDEs. 
This result was generalized in m in d-dimension and by considering the class of one-sided Lipschitz 
drift coefficients. The proof of results in mm are based on a Yamada-Watanabe approximation 
technique (see [2T]). In the jump-diffusion case, the work [T9] shows that when the coefficients are 
non-Lipschitz, the Euler-Maruyama approximation has a strong convergence. The concept used 
to prove their result is to the Yamada-Watanabe approach and is based on a generalized Gronwall 
inequality. 

Let us also mention the work [13] , where the author prove that if L is a Wiener process and the 
coefficients are /3-Holder continuous, then the Euler-Maruyama approximation converges weakly 
to the unique weak solution of the corresponding SDE with rate /3/2. This result was extended to 
the case of nondegenerate SDEs driven by Levy processes in [14]. More specifically, the authors 
study the dependence of the rate on the regularity of the coefficient and the driving noise. The 
backward Kolmogorov equation plays a crucial role in their argument. 

In this paper, we first to study the strong rate of convergence of the Euler-Maruyama scheme 
(|2|) when L is a Wiener process and the drift coefficient b is singular, that is b is bounded and 
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/3-Holder continuous in space variable. Our method to study the strong rate of convergence differs 
from the existing ones. We develop a method based on the regularity of the solution to the 
Kolmogorov equation associated to the SDE dT|) . More precisely, using Ito’s formula, we write the 
drift part in terms of the solution to the Kolmogorov equation. Applying the estimates of the 
derivatives of the solution to the Kolmogorov equation, we are able to obtain the convergence rate 
(see Theorem l2.11l) . 

Second, we examine the Euler-Maruyama scheme m when the driving noise is a truncated 
symmetric cc-stable process with a £ (1,2). The method is also based on the regularity of the 
solution to the Kolmogorov equation associated to the SDE ([!]). The authors are not aware of any 
other work where the regularity of the Kolmogorov equation is used to study strong convergence of 
Euler-Maruyama approximation with irregular coefficient. Traditionally, the Kolmogorov equation 
based method is used to study the weak convergence rate for the Euler-Maruyama approximation 
(see for example [IBS and references therein) which is for example very important in financial 
applications. 

This paper is divided as follows: In Section [2j we introduce some notations and preliminary 
results on existence and regularity of the Kolmogorov backward equation associated to the SDE 
CD- The main results of this paper are also given in this section. Section [3] is devoted to the proof 
of the main results. 


2 Preliminaries 


2.1 Notations 


Let jB(Mq) be the Borel-cr-algebra on Rj), with Rg := R d \ {0} and set V = D = ..., gf^)*, 

D 2 = ( dxix- an d A = Eti H ere * is the transpose of a vector or matrix. In the 

following we introduce some space of function: 


• Cb(R d ;R fc ), d, k £ N denotes the space of bounded continuous functions from R d to R fe . In 
particular, if k = 1, we say Ct,(R d ;R) = Ct,(R d ). For bounded measurable function /, the 
supremum norm of / is defined by ||/||oo := sup^ggd |/(x)|. 

• C'£°(]R d ) denotes the space of all infinitely differentiable R-valued functions with compact 
support contained in R. 

• (R d ; R fc ), f3 € (0,1) denotes the set of all functions from R d to which are bounded and 
/3-Holder continuous functions. Hence if f £ Cf (R d ;R fc ), then 


sup 

x,yEM. d ,x^y 


1/0*0 - f(y) I 

\x-y\P 


< oo. 


• Cj ,,3 (R d ), i = 1,2 and /3 £ (0,1) denotes the space of /-times differentiable functions / : 
R d —> R with D i f £ C^(R d ;R®^) for any 1 < i < i. A function / : R d —> R d belongs to 
(7^(R d ;R d ) if each it components fj is in C^( R d ) for j = 1,... ,d. 
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• Let F be a class of functions and [a, b] be closed interval. C{[a, b\, F) denotes the space of all 
functions / : [a, b] x R d —> R fc such that f{t, •) £ F for any t £ [a, b\. 

• C 1 ([a, 6], F) denotes the space of all functions / such that / £ C([a, 6], F) and |£(i, ■) exists, 
is continuous and is in F. 


• For a < b, we write Cf([a, 6]) for C([a, 6]; Cf (R d ; R d )) and define the norm || • || c /3([ ob p on 
Cf ([a, b]) by 


'C§([a,b]) 


sup \f(v,x)\ 


sup 

v€[a,b] ,x^y 


\f(v,x) - f(v,y )I 
\x - y\P 


2.2 Levy processes: definition and basic properties 


In this section, we recall the definition and basic properties of a Levy process. We also give the 
definition of truncated symmetric a-stable process. For more information on Levy processes, we 
refer the reader to pom Emm and references therein. 

Definition 2.1. A stochastic process L = ( L t )o<t<T on R d is called a Levy process on probability 
space (f2, T, P) if the following conditions are satisfied: 


(i) For any choice of n £ N and 0 < t\ < ■ ■ ■ < t n , the random variables L to ,L tl — L to ,. .., L tn — 
Lt n _ x are independent. 

(ii) L 0 = 0, a.s.. 

(iii) For any t, s > 0, the distribution of L s+t — L s does not depend on s. 

(iv) For any e > 0 and t > 0, lim s _,. t P(|L S — L t \ > e) = 0. 

(v) There exists tto £ T with P(flo) = 1 such that, for every u £ floj L t (w) is right-continuous 
in t > 0 and has left limits in t > 0, i.e. L t (uj) is cadlag function on [0,T]. 


Next, we define a Poisson random measure associated to a Levy process L. Let A L t := L t —L t - 
be the jump size of L at time t. We define a Poisson random measure for L on S([0, oo)) x B(Rg) 
by 


N(t,F) := J2 1f(AL s ), FGB (<). 

0<s<t 


Then the Levy measure v is defined by v{F) = E[AT(1, F)) for F £ Z?(Rg). 

By the Levy-Ito decomposition (see for example cm Theorem 19.2]), a Levy process L admits 
the following integral representation: 


Li — bt + oWt T 


0 J\z\<l 


zN(ds , dz) + 


o J\z\>\ 


zN(ds , dz), 


where b £ R d , cr £ R d x R d , W = {W t )o<t<T is a d-dimensional Wiener process and N defined by 

iV(ds,dz) = 7V(ds,dz) — ^(dz)ds, 


is the compensated Poisson random measure of L. 
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Remark 2.2. Let L = ( L t )o<t<i be a d-dimensional Levy process with E[|L t | p ] < oo for some 
p > 2. Then it follows from EB Theorem 66] that, there exists C v > 0 such that for all t £ [0,1], 
we have 


E 


sup \L S \ P 

,0<s<t 


< C p t. 


A d-dimensional (d > 2) truncated symmetric a-stable process L with a 6 (1,2) is a Levy 
process such that 

E[e i < f " t *>] = e - *^ 0 , (eR d . 


Here the function ip is the characteristic exponent of L and is given by 


'!’(€) = A{d,-a) [ {1 - cos((£,z))}i/(dz), 

J\z \<1 


with 

A(d, —a ) := 


and the Levy measure given by 


a2“" 1 r((d + a)/2) 

7r d / 2 T(l - a/2) 


i/(dz) 


i(kl < i) 


I d+a 


dz. 


It follows from the Levy-Ito decomposition, a truncated symmetric a-stable process L admits 
the following integral representation: 


L t 



( 3 ) 


Note that for any n £ N, \z\ n v(dz) = 0. Then using [TJ Theorem 2.5.2], the n-th moment of 
L t is finite for any t > 0. In particular, the second moment is finite. This means that a truncated 
symmetric a-stable process L is a square integrable Levy process. By the integral representation 
([3| and the Ito isometry, we conclude 


E[|L t | 2 ] = f [ \z\ 2 v(dz)ds = t [ 
Jo d|z|<l J\i 


1 


N<i 


\d-\-a — 2 


zwd y = ct, 


for some constant c. 


Finally the infinitesimal generator £ of a truncated symmetric a-stable L is given by 
Cf{x):= [ (f(x + z) - f{x) - {z,Df{x)))v(dz), f&C™{R d ). 


2.3 Results on existence of solutions to the Kolmogorov equations 

In this section, we present some auxiliary results on existence and regularity of solutions to the 
Kolmogorov equations. Further, we derive some explicit expressions of the constants in the upper 
bounds of the derivative of these solutions. These results will play a crucial role in the proof of 
our main theorems. 
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2.3.1 Case of Wiener process 

The following result corresponds to [U Theorem 2.8] provides regularity of solution to the Kol¬ 
mogorov equation associated to the SDE ([!]) when the driving process is a Wiener process. 

Theorem 2.3. Assume that L = W be a d-dimensional Wiener process. Let t\ £ [0,T]. For 
all ip £ C([0, ii]; C/ (R d ; R d )), there exists at least one solution u to the backward Kolmogorov 
equation 

du 1 

+ Vu ■ b + -Am = 011 [0, ti] x u(ti,x) = 0 (4) 

of class 


u £ C([0, <r]; C 6 2l/3 '(R d ; R d )) n C' 1 ([0, fr]; cf (M d ; R d )) 


for all 0' £ (0,(8) with 

H £> 2 “Hcrf([0,t 1 ]) - ^'Ikllcfao.ti]) 

and 


l|V«llcf ([ o,t l]) ^ C-(*i)ll^llcf([o,t l] ) with = °- 

In the following Lemma we make precise the function C(ti) derived in Theorem 12.31 

Lemma 2.4. The constant C(ti) defined in Theorem 12.31 is given by C(ti) = Cot / for some 

1 /2 

constant Co and £ (0,T] with \\b\\ c f<^ 0 T ^Cot{ < 1/4. 


Proof. Let us first prove that if v is the solution to the heat equation 

dv 1 

= - Av + tp on [0, ti] x R d , v(0, x) = 0, 
with ip £ C([0, ti]; Cfr (R d ; R d )), then it holds that 

ll Vw llcf([o,ti]) - C(*)IMIcf ([0,ti])’ 

It is shown in 2) Theorem 2.3] that 


j(t,x)= [ E[ip(s,x+ W t -s)]ds,t £ [0,t! 

Jo 


is a solution to the equation J5]) and 


-^—E[(p{s, x + W t -s)} = -j^E[cp(s, x + W t -s)Wi_ s }, 


( 5 ) 


( 6 ) 
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where W{ is the i-th coordinate of W t - Since the function \x\ exp(—|a;| 2 ) is bounded, there exists 
Co > 0 such that 

|Vu(£,x)| < C 0 sup \ip(s,x)\ - ds < C 0 t\ /2 \\ip\\ c p, [ ]v (7) 

se[o,ii],xeR d Jo y/t-s 6 u 1J ' 


where Co = 2Co- Similarly, we get 
\X7v(t,x) - \7v(t,y )| 


'-y\ p 


< Co sup 

se[0,ti],x^y 


|<p(s,x) - <p(s,y )| 


y|= d * s c » ( l / 2 ll* , llo< i ([o,t l] )' 


( 8 ) 


Combining © and (| 8 |). we conclude that © holds. 

Next, set = 0 and for n £ N U {0}, let be a solution to 

o„.( n + 1 ) i 

+ -A U (” +1 ) = -(&• V)u w - cp on [ 0 ,ti] x R d , «( n+ 1 )(t 1 ,s) = 0 . 
Define v^ n \ n £ N U {0} by v^ := (b ■ V)u^ + tp. From ©, we have 

l|V« , ” + 1 , ll o » ( | 0 ,„ l) <C(t 1 )l|V»‘“>|| cf(I0 ,, l|) 
for any n £ N U {0}. Choose t\ with || 6 || c 3 q 0 t ^C(C) < 1/4, we have 


l|V« (n) llc?([0, t i]) — IMIcf([o,ti]) + 2 || 6 || cf ([ 0 it i])||V^)|| c ,( [ 0 it i]) 

— IMIcf ([0,t]) + 2 ll 6 llcf([0,T]) C '( tl )l|Vw (n ([0,ti]) 

n 1 

< « * < X] 2» H^llcfao.ti]) < 2 IMIcf ([o,td)- 

i=l 


Hence, it holds that 


l|VM ( ' !+1) || c ^([ 0 , tl ]) < 2C(ti)||(/?|| C f!([ 0 ,ii])■ 

On the other hand, we know from [4j Theorem 2.8] that there exists a subsequence (u^ nk ^)keN of 
(u ^) n &n that converges uniformly in (£, x) £ [ 0 , £i] x R d to some u such that u is a solution to the 
backward Kolmogorov equation © and V'«’■”'=) converges to Vu in Cf([0, £i]). Hence the result 
follows. □ 


Similar arguments as in the proof of Theorem 12.31 and Lemma l2~4l lead to the following corol¬ 
lary: 

Corollary 2.5. Let 0 < £i < £2 < T with ||^|lcf([o t])Co ■ (£2 — ti) 1 ^ 2 < 1/4. For all ip £ 
C([£i, £2]; C p (K d ; R d )), there exists at least one solution u to the backward Kolmogorov equation 

du 1 

■ + Vu • b + - A u = —p on [£ 1 , £ 2 ] x R d , «(£ 2 , x) = 0 
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of class 


u £ C([ tl M:C 2 /\M. d - R d )) n C 1 (fr, i 2 ]; cf(R d ;R d )) 
for all /?' £ (0, (3) with 

H £> 2 ' u llcf'([ti,t 2 ]) - M Mc^m) 

for some constant M and 

l|Vu|lcf([ t i, t2 ]) - C o(*l _ *2) 1/ " ll^llof ([ti,i 2 ]) 

for some constant Co- 

The following result will also be needed in the proof of the main results. 

Lemma 2.6. Let T > 0. For any e £ (0,1), there exist to £ N and {Tj)j -such that 
0 = Tq < Tj < Tj + 1 < • • • < T m = T and for any j = 0, ..., to — 1, 

IMIcfao.TpC'o • ( T i+1 - T i) 1/a < e and ll & llcf([0,T]) C 'o • (T j+1 - Tj) 1 / 2 < i. 


'C£([0,T]) 

Proof. For any e £ (0,1), we define 


<5 := 


k ^oll^llcf ([0,T]) ([0,T]) / 


Since 6 > 0, there exists to £ N such that (to — 1)<5 < T < mS. Dehne To := 0, T m := T and 
Tj := Tj _i + (5 = jS for j = 1,..., m — 1. Then for any j = 0, ..., m, Tj — T,_i < 5, we have 

IMIcf ([o,t])^o • (Tj+ 1 — Tj) ! < e and 11 ^ 11 ([o, x 1 ] ) ~ 7 < 

This concludes the proof. □ 


2.3.2 Case of truncated a-stable process 

The following result is due to no Theorem 17] and it provides existence and regularity of the 
Kolmogorov equation in a given space when the driven noise is a truncated symmetric a-stable 
process. 

Theorem 2.7. Assume that L = ( L t )o<t<T is a d-dimensional truncated-a-stable process for 
a £ (1,2) and d> 2. Let ti £ [0,T]. Suppose that ip £ C([0, ti]; Cf (R d ; R d )) for (3 £ (0,1) with 
a + (3 > 2. Then, there exists a u £ C([0, ti], C 2 (R d ; R d )) D C 1 ([0, ti], C^R 6 *; R d )) satisfying the 
backward Kolmogorov equation 

du 

— + Vu • b + Cu = —ip on [0, t\) x R d , u(t\,x) = 0, 


( 9 ) 




with 


ll Vu llcf([ 0 )tl ]) ^ C'«(*i)|| < ^Hcf([o, tl ]) with = ° 


and 

\\D u||oo < MM c e {[o ti]) , 

for some positive constant M. 

The function C a (t\) is made more precise in the following lemma: 

Lemma 2.8. The function C a (ti) dehned in Theorem 2.8 is given by C a (ti) = Cot\ for some 
constant Co and t\ S (0,T] with IHIc^qo < 1/4. 

Proof. Let us first prove that the solution v to the equation 


du 

— = Cv + ip on [0, <i] x R , u(0, x) = 0, 
for p £ C([0, ti]; C^(R d ; R d )) satisfies 

ll Vw llcf([o,ti]) — < ^a( il )Hi 9 llcf([o.ti])- 
It follows from [6j Theorem 16] that 

v(t,x)= f E[<p{s,x + L t - S )]ds,t S [0,fi] 

Jo 

is a solution to the equation (flOl) and (see (19) of 0) 

o /» O 

—E[tp(s,x + L t )\ = - / p(s,t 1/a u +x^^—— (p t {t 1/a u))di 

OXi J R d OXi 


( 10 ) 


( 11 ) 


Here Pt{x) is the density function of L t which satisfies 


IR d 




d u < C' 0 t~ d / a 


for some constant Cq , (see, Hj , page 5338). Hence there exists a constant Co > 0 such that 


\Vv{t,x)\<C 0 sup \p(s,x)\ I ——l-^ds^ C 0 t\ 1/ “|b|| c P ([otl]) , (12) 

se[o,ii],xeR d Jo \t~ s ) b ’ 


where Co := a/(a — l)Co- Similarly 
\Vv(t,x)-Vv(t,y)\ ~ \<p(s,x)-p{s,y)\ rt 

— \ L/q Slip 


-y\P 


«6[0,ti ],x^y I® V\^ Jo (1 ®) 


J 0 (t - s)i/“ ds - c ° tl 1 ([0,ti])- 

( 13 ) 
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Combining (fl2l) and m, we conclude that m holds. 

Next, set u^ = 0 and for n £ N U {0}, u( n+l ' > be a solution to 


du ( n+1 ) 
dt 


+ £u (n+1) = -(&• V)u {n) -<p on [0,fi] x R d , u {n+1 \t i,x) = 0. 


Arguing as in Lemma l2~4l we get by choosing t\ with ||&|| c ^ 0 T ^C a {t\) < 1/4 

ll W ’ l+1) |lcf([0, t i]) - 2 ^ 7 “( il )ll^llc/([0,ti])' 

Furthermore, we know from [6] Theorem 17] that the sequence (u^ n ^) n6 N which converges uniformly 
in (t,x) £ [0,ii] x to some u such that u is a solution to the backward Kolmogorov equation 
© and Vu^ converges Vu in Cf([0, ti]). Hence the result follows. □ 


The following result can be derived using similar arguments as in the proof of Theorem 12.71 
and Lemma [2~8l 

Corollary 2.9. Let 0 < t\ < t 2 < T with ||^|| C 7' 9 ([o t])^° ■ (t 2 — <i) 1_1 / a < 1/4. For all <p £ 

C([ti, t 2 \; Cfr (M d ; M d )) for /3 £ (0,1) with a + f) > 2, there exists a u £ C([ti, i 2 ], (R d ; R d )) fl 

C 1 ([ti, i 2 ], Cfc (R d ; R d )) satisfying the backward Kolmogorov equation 

du 

— + Vu ■ b + Cu = —p on [fi,f 2 ] x u(t 2 , x) = 0, 


with 


and 


||Vu|| 


c£([ti,t a ]) 


<Co(t 2 -ti) 1 - 1/ “||^|| 


c£([t iM) 


\\D 2 u\\oo < MM 




for some constant M. 

We also have 

Lemma 2.10. Let T > 0. For any e £ (0,1), there exist m £ N and (Tj)j—o m such that 
0 = T 0 < Tj < Tj + i < • • • < T m = T and for any j = 0,..., m — 1, 

IMIc£([o,T])Co • (T j+ 1 - Tj ) 1 - 1 /“ < e and ||6|| cf([0iT]) Co • (T i+1 - Tj) 1 ^ < ± 

2.4 Main theorems 

In this section, we state the main theorems of this paper. We obtain results on the rates of the 
Euler-Maruyama approximation in L p -sup norm for p > 1 if L is a Wiener process or a truncated 
symmetric a-stable process. 
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Theorem 2.11. Let L = W be a d-dimensional Wiener process. Assume that the drift coefficient 
b is bounded /3-Holder continuous with /3 £ (0,1) in space and 77 -Holder continuous in time with 
77 £ [1/2,1], i.e., there exists K > 0 such that for any x, y £ M. d and £, s £ [0, T], 

| b(t, x ) — b{t , y)| < A"|x — y\& and | b(t, x) — 6 (s, x)| < K\t — sj''. 

Then for any p > 1, there exists a positive constant C depending on K, M, T, d,p, xo, /3, 77 and 
ll^llcf ([0 t]) SUC ^ that 


E 

sup 

x t - x t (n) 

P~ 


0 <t<T 


_ 


C 

— n pp/ 2 ' 


Remark 2.12. Under conditions of Theorem l2.111 the pathwise uniqueness holds for the SDE CU 
(see , [41 Corollary 2.3]). Then [51 Theorem D] guaranties that the Euler-Maruyama approximation 
@ converges to the corresponding SDE in L 2 -sup norm. Hence, Theorem 12.111 generalizes [5} 
Theorem D] in two directions: first, it gives an L p -sup convergence for any p > 1 and second, it 
gives the rate of convergence which was not given previously. 

Theorem 2.13. Assume that L is a d-dimensional truncated symmetric a-stable process with 
a £ (1, 2) and d > 2. Suppose that SDE ([I]) has a unique strong solution. Moreover, assume that 
the drift coefficient b is bounded, /3-Holder continuous in space with /3 £ (0,1) and a + (3 > 2, and 
77 -Holder continuous in time with 77 £ [1/2,1], Then, for any p > 1, there exists a positive constant 
C depending on K, M, T, d, p, xq, a , /3 ,77 and ||d|| c ./ 3 q 0 such that 


E 


sup 

x t - x[ n) 

V 

0 <£<T 


. 


G 


-iP0/2 


if p/3 > 2 , 

if p > 2, 1 < p/3 < 2 or p £ [1, 2). 


Remark 2.14. Let L be a d-dimensional truncated symmetric a-stable process with a £ (1,2). 
Choose the drift coefficient b of the form 


m 

*= 1 

where /» are continuous functions and bi are as in Theorem 12.71 for i = 1 ,...,m. Then the SDE 
dT]) has a unique strong solution on small interval [0, T] (see [HI Theorem 23 and Remark 24]). In 
Remark 1 2.2 1 we notice that for all L p -integrable Levy processes for some p > 2, the p-th moment 
is always bounded by C p t. Therefore, the rate of L p -convergence with p > 2 and p/3 > 2 coincides. 


3 Proof of main theorems 

This section is devoted to the proof of the main results. The constants C\, C 2 , C 3 and C are assumed 
to be positive and independent of n. Unless explicitly stated otherwise, the constants C\,C 2 ,Cz 
and C depend only on K, M,T, d,p,Xo,a, /3, A, 77 and || 6 || c / 3 q 0 T j). Moreover, the constant C may 
change from line to line. 
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We define 


Xt := {Xl ..., X?)* and X t (n) := (X^ n \..., x\ d ' n) )*. 


The following estimation is standard. For the convenience of the reader, we will give a proof. 
Lemma 3.1. Suppose that the drift coefficient b is bounded and measurable. 

(i) If L = W is a d-dimensional Wiener process, then for any p > 0, there exists C > 0 such 
that for any t £ [0, T], 


E 


xi n) - x 


(n) 

Vn(t) 


p 1 


< 


c 

T,P/ 2 ' 


(ii) If L is a d-dimensional truncated symmetric a-stable process with a £ (1, 2) and d > 2, then 
for any p > 0, there exists C > 0 such that for any t £ [0, T], 


E 


xi n) - X 


(n) 

Vn(t) 


< 


c 


n P/2 


if P > 2, 
if p < 2. 


Proof. From the definition of the Euler-Maruyama approximation, it holds from boundedness of b 
that 




r ( n ) 


T\\b\\ c 


+ | L t -L 




Suppose L = W. Then for any p > 1, 

pi 2P- 1 TP||6||g 0 


E 


v( n ) V 

X t X r,n(t) 


"(n) 


< 


nP 


2 P-1 E [\W t 


W, 


J — 


< 


c 

n p / 2 ’ 


Hence the statement is true in the case of Wiener process. 

Suppose that L is a truncated symmetric a-stable process. It is enough to prove the statement 
for p> 2. Since E[|L t | p ] < oo it follows from [17] Theorem 66] that there exists C such that 

E[\L t -L, n{t) n=E[\L t . Vn{t) \ p ]<^. 

Hence the result follows when L is a truncated symmetric a-stable process. □ 


3.1 Proof of Theorem 12.111 

For a given e £ (0,1), we consider the partition of closed interval [0,T] which is 

considered in Lemma 12.61 For i = 1,..., d and j = 1,... ,m, Corollary 12.51 implies that there exists 
at least one solution Uij to the backward Kolmogorov equation: 

du 1 

+ Xuij ■ b + - A = -bi on [Tj-i,Tj\ x K d , u it j{Tj,x) = 0 
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and Uij satisfies, 


ll^^dllcf[ t 3 -_i,t 3 -] — Co • {Tj Tj_ i) ! ||6jllc f ^([Tj_i,Tj]) — Co • (Tj Tj-i) ^ IHIcf([o,T]) — e - 

For any i e by Ito’s formula, we have 


'U'iij (t: Xt') — Uij Xxj-x ) “t - 


<9lt 


^ (s,X s )ds + J S7uij(s,X s )dX s + - 


Auij(s,X s )ds 


= Ui i j(T j - 1 ,X Tj _ 1 ) - f bi(s,X s )ds+ [ Xu i} j(s,X s )dW s . 
JTj-i -M-, 


Hence we have 

r t 


[ bi(s,X s )ds = Uij(T j _ 1 ,X Tj _ 1 ) - Uij(t,X t ) + [ Vu id (s,X s )dW s . 

JTj-i JTj_i 

In the same way, we have 

f bi(s,X^)ds = UijiTj^X^ J - UiJ (t,X^) + [ Wu itj (s,X^)dW a 
J Tj -1 3 JTj-i 

+ f Vutj(s,X^) ■ (bi(Vn(s),x£\ s) ) - bi(8,xWj) da. 

JTj -1 

It follows from (fill) and (fT5l) that for any i = 1 ,d, 


\ri xr(n,i) _ v{ n 

A t ~ A t — A T,-_ 1 — A T. 


C + r (6,(5, X s ) - bi(r) n (s), X$ 8) )) ds 
JTj -i V 


= X^ - X, 


(«.*) 


3-1 lj-1 


+ (ujj (Tj -!, Xy J _ 1 ) - ^(Tj-i.xWj) - ( Ui)i (t,X t ) - ^(t,X t (n) )) 
+ (v«ij(a,X.) - Vujj(s, xW)) dW s 

+ f VuijfaxW) • ( 6 ( 5 , X 8 W) - 6 (^( 5 ), X^ (8) )) ds 
+ f (&,(s, Xj">) - 6 *^( 5 ),X^ w )) ds. 

JTj-j 

Since UVujjH^Qj,. i T .^ < e, by the mean-value theorem, we have 


*t* - 


(n,i) 


< 


vi v{ n J) 

A 71 . _ Am 

ij-1 1j—1 


l i,j(Tj- 1 , Xxj-j) Uij(Tj- 1 , Xjx^) + Ujj(t, X t ) Uij (t, Xj 1 ) 




T,- : 


(v Ui j(s,X.) - Vuijfs.xW)) d w s 


(14) 


(15) 
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+ f VuijfaxW) b(s,X^)-b( Vn (s),X^\ s) ) 

•' T i -1 

■’Tj -1 

< (1 + e) 


ds 


ds 


It,,, - X, (n) 


i - 1 Tj_i 

t 


x t - x f (n} 


fVu itj (s,X s ) - Vuij(8,XW)) dW s 

r j -1 V 

+ (1 +e)K £ ^ + ^+^)K(t-T j _ 1 )(F.y . 


For p > 2, using Jensen’s and Holder inequalities, we have 


p/2 


Xt-X t 


(n) 


= Ek>‘- X , 


(n,i) 


,i=1 


< dP/ 2 - 1 Y, I Xt - X, 


(n,i) | 


i=l 


x t -x t 


< d p/2 5 p ~\ 1 + e) p X Tj _ 1 - X^Y + d p/2 5 p ~ 1 e p 

+ d P/ 2 -i 5 p -iJ2 I hu ilj (s,X a )-Vu i ^(a,X^)\dW l 

i =1 dTj -1 V 

rt 

+ d p/2 5 p -\l + e) p K p {t - Tj-if - 1 / Xi n) - X^ {s) 

JTj-1 


(n) 


P0 


ds 


+ d p/2 5 p -\ 1 + £) p X p (f - 


Since e > 0 is arbitrary, it may be chosen such that c(p, d, e) := d p / 2 1 5 p 1 £ p < 1. Then we have 


J-l 

PV 


x t -xy 


p d p ! 2 SP-^l +£) p 
(1 -c(p,d,s)) 

d p/2-l 5 p-l d 

+ (l-c(p,d,£)) ^ 




'Ti-, 


(vu id (s,x s ) - VByfs,!'" 1 )) dtr s 


+ £ y K P( t _ T J _i ) p - 1 
(1 - c(p, d,£)) 

d p / 2 5 p ”i ( 1 + £) p X p (t-T j _ 1 ) p ^T 
(1 - c(p, d, e)) 


Tj-x 

PV 


xl n) - x 


Pn(s 


p/3 


ds 


(16) 


Taking the supremum and then expectation on both sides of (|16H , we have from Burkholder-Davis- 
Gundy’s inequality and Jensen’s inequality that 


E 


sup 


X - 

x\ u 


< d p/2 5 p ~ 1 (l+£) P E 
(1 -c(p,d,£)) 


y _ y( n ) 
A T j -1 


rf p/2-l 5 p-l C ( P;d ) r §-l d rt 

(1-C(p,d,£)) 


Vu j]i (s,X s )-Vi li|j (s,xW) 


ds 
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dP/ 2 5 p -i(l. + e)PK p TP~ 1 
(1 - c(p,d,e)) 
d p/2 5 p-i(i + e y(KT)P f T 


E 


Tj -1 

VP 


x< n > - x ( 


Vn(s) 


P$ 


ds 


(1 -c(p,d,e)) \n 

where C(p,d) is the constant in Burkholder-Davis-Gundy’s inequality. From Lemma 13.11 (i), we 
have 


E 




p 

rfp/ 2 5P - 1 £P 


p-1 

sup 

Tj — i <n<t 



< E 

(1 - c(p,d,e)) 

u 

l 

- 


dP/ 2 5 p ~ 1 C{p,d)T^- 1 £P 
(1 - c(p, d,e)) 

d p / 2 5P-\ 1 + e)P(KT)P C + Tpp 
(1 — c(p,d,e)) riPP/ 2 


[ E 

sup 

X - x(”) 

p 

M'-1 

Tj — i<.u<.s 




ds 


CiE 

£ 

1 

1 

p- 

+ C 2 [ E 

sup 

x - 

P 


- j. 


JTj_i 

Tj—i <u<s 




Next, we prove by induction that for each j = 1,..., m, 

As 


E 

sup 

X - xW 

p 


Tj_i <n<t 




< 


j.P/3/2 


, t e (Tj-i,Tj], 


ds 


(17) 


C 3 


nP/3/2 ' 


where A\ := C 3 e C2T and Aj := {C\Aj_\ + C 3 )e C2T for j = 2,_ ,m. If j = 1, since To = 0, for 

any t £ (0,Ti], we have 


E 

sup 

x - x^") 

V 

<C 2 f E 

sup 

x - x 

p" 


_0<u<t 



Jo 

_0<u<s 


_ 


ds + 


C 3 

n P/3/2 ' 


By Gronwall’s inequality, we have 

E 


sup 

x — x(") 

p" 

_0<n<£ 


. 


< 


C 3 e c * T 


n P/3/2 

Assume that holds for j = 1 ,...,i — 1 with 2 < i < m. Then, for any t £ (Tj_i, Tj], we have 

ds + - 


E 

sup 

X - X*t) 

p 

< CiE 

Xti-i - 

P- 

+ C 2 [ E 

sup 

x - x(") 

P 


Tj_i <n<t 






JTi- 1 

Tj_i <u<s 




G 3 


7 p/3/2 


<C 2 E 

JTi -1 


sup 

X - x( n ) 

P 

Ti—i <u<s 




ds + 


CiA,.! + C 3 


7 p/3/2 


Using once more Gronwall’s inequality, we have 


E 

sup 

X - X*-"^ 

P 


Tj_i <n<t 




< 


+ C 3 )e 2 A. 


nPl 3/2 


tP/3/2 ' 
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Hence dD holds for any j = 1,... ,m and we have 


E 

sup 

x s - xW 

p" 

m 

<E E 

sup 

x - W”) 

P 


-0 <S<T 



3=1 

Tj — i <u<Tj 




This concludes the statement of Theorem 12.111 



3=1 


□ 


3.2 Proof of Theorem 12.131 


Let us first remark that for p £ [1,2), the L p - norm is bounded by L 2 -norm. Hence it is sufficient to 
prove the statement for p> 2. As Theorem 12 .Ill let e £ (0,1) be given and consider the partition 
(Tj)j- o,...,m of closed interval [0, T] defined in Lemma 12.101 For i = 1,..., d and j = 1,..., to, 
Corollary 12.91 implies that there exists at least one solution Uij to the backward Kolmogorov 
equation: 

du ' ’ 

+ Vu it j ■ b + Cuij = -bi on [Tj- i, Tj] x R d , u it j(Tj,x) = 0 

and Uij satisfies, 


ll^billcf ([Tj_i,Tj]) — Co(Tj 


rn \ 1 — l/a: 

1 j-1 ) 


INIcfa Tj-uTi]) 


< Co (Tj 


rri \1 —I/O! 
1) 


cfao.T]) C e 


and 

\\D Mijlloo ^ ([0,T])' 

For any t £ [Tj_i,Tj\, using Ito’s formula and Kolmogorov equation, we have 

/ bi(s, A s )ds — 'U'itj (lj_l, X-Tj ) ,j (t, Xf] 

JT }-1 

ft j. 

{Uij(s, A s _ + z) - Uij(s,X s _)}N(ds,dz) 


,T j -1 J\A<d 


and 


t h(s,X^)ds = UijiTj-uX™) - Ui tj (t,xl n) ) 

•'Tj -1 

+ [ [ {Ui,j(s, X ( ^ + z) - u itj {s , X^)}N{ds, dz) 

J Tj-l J\z\<1 

+ f VuijfaxW) ■ (b( Vn (s),X^ (s) ) - 6 (a| n) )) da. 

•'Tj -1 


(18) 


(19) 


( 20 ) 


Combining (19| and (1201) and using similar arguments as in the case of Wiener process we have, 


x t -x t 


(n) 


P <p/ 2 b p - l {\ +e) f 
(1 -c(p,d,e)) 


v y( n ) 

A T-i a t,-_ i 
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+ 


d p/2- l 5 p— 1 d 

(1 -c(p,d,e)) j T ._ 1 j | z[5 i 

dP/ 2 5 p - l {l + e) p Kv{t-T j _ 1 y r t 

(1 -c{p,d,e)) J T ._ X 

d p/2 5 p-i( 1 + e yKP{t - Tj_i) p pv 

(1 - c(p, d,e)) 


F i j(s,2)iV(ds ! d0) 


X (n) _ X (n) 

S ^n(s) 


d s 


( 21 ) 


where 


Hij(s, z) := {u i:j (s,X s _ + z) - iti,j(s,X s _)} - |iijj(s, + z) - u itj (s, X s ( " } )}. 
Using twice the mean value theorem, we have 


d r i 

_ \ ' / ( yk \r(n,k)^ 


HiAs,*) = ±: I (X*--Xr- K> ) 


du. 


dx k 


?(»• 


Xt ] + « + r(X s _ - x£ } )) - (s,X^ + r(X s _ - x£ } )) !> dr 


j\ej\r ±(xl - (x« + r(*,_ - *£?) + to) . 


From Kunita’s inequality (see [2 Theorem 4.4.23]) and (fTHll . we have 


E 


sup 


r / 

J Tj-l J\z\<l 


z)N(ds, d z) 



/ /■* r \ p/2 


/•t n 

< CE 

/ / | J ff I ,j(s,z)| 2 I /(dz)ds 

+ CE 

/ / |i?i,j(s, z)| p i/(dz)ds 


V/t,--! J|*|<1 / 


JTj_i J|z|<1 



/ rt r 


2 \ P/2 1 


rt r 


< CE 

V'U-i -'bl<i 

x s - X W 

|z| 2 z/(d 2 )ds 

+ UE 

JTj -1 i|z|<l 

x s - x(") 


[ E 

x s - X W 

p- 

•Tj -1 




ds. 


/(dz)dfi 


( 22 ) 


Taking the supremum and then the expectation on both sides of (1211) and using (1221) and Lemma 
13.11 (ii), we have 


E 


sup 

Tj_ i <u<t 


X - X (n ) 


< CiE 

£ 

1 

1 

>• ^ 

P- 

+ C 2 [ E 

sup 

x - x( n ) 

p 


- j 


•'U-i 

Tj — i <n<s 




ds 


+ C / E 

JTj-\ 


x) n) - x 


»?»(«) 


p/3' 


c 

ds H- 

u pv 
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< C lE 


if p/3 > 2, 
if 1 < p/3 < 2. 









1 

< CiE 

$ 

1 

1 

& 

P- 

+ C 2 f E 

sup 

i 

4 

P 

ds + < 


j 



Tj — i <u<s 





C 3 



n pp/2 


Arguing as in the proof of (fI71) . we have that for t G (Tj_i,Tj], 


r 




sup 

X - A (n) 

P 

<\ 

Tj_i <n<t 









Aj 



n p P / 2 


if p/3 > 2, 
if 1 < P/3 < 2, 


where Ai := C^e c,2T and Aj := {C\Aj-i+C^)e c,2T for j = 2 ,... , to. This concludes the statement 
of Theorem l2.13l □ 
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